plot(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=2)
grid(equilogs = T, col="darkgrey")
abline(v=prop, lty=5, lwd=1.5)
Plot_ErrorBars(x = mutryk[,1],y = mulokc[,1],x.err = sdtryk[,1],y.err = sdlokc[,1],xbar = T,ybar = T)
Plot_ErrorBars(x = mutryk[,2],y = mulokc[,2],x.err = sdtryk[,2],y.err = sdlokc[,2],xbar = T,ybar = T)
points(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,40), ylim=c(10,40), pch=16, col=2)
points(mutryk[,2],mulokc[,2], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=3)
fit8 <- lm(mulokc[,1]~mutryk[,1])
fit16 <- lm(mulokc[,2]~mutryk[,2])
summary(fit)
abline(fit8, col=2, lty=2)
abline(fit16, col=3, lty=2)
coef8 <- coefficients(fit8)
coef16 <- coefficients(fit16)
sig <- 3
text(48,30,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=.8)
text(48,27,labels = paste("16 mm: ",signif(coef16[2],sig),"x+",signif(coef16[1],sig), sep=" "),col=3, cex=.8)
legend("topleft",legend = c("Proportionerede styrker","Max. tilslag: 8 mm", "Max. tilslag: 16 mm"), lty=c(5,0,0), pch=c(NA,16,16), lwd=1.5, col=c(1,2,3))
# ----------------
par(pty="s",xpd=TRUE)
plot(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=2)
grid(equilogs = T, col="darkgrey")
abline(v=prop, lty=5, lwd=1.5)
Plot_ErrorBars(x = mutryk[,1],y = mulokc[,1],x.err = sdtryk[,1],y.err = sdlokc[,1],xbar = T,ybar = T)
Plot_ErrorBars(x = mutryk[,2],y = mulokc[,2],x.err = sdtryk[,2],y.err = sdlokc[,2],xbar = T,ybar = T)
points(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,40), ylim=c(10,40), pch=16, col=2)
points(mutryk[,2],mulokc[,2], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=3)
fit8 <- lm(mulokc[,1]~mutryk[,1])
fit16 <- lm(mulokc[,2]~mutryk[,2])
summary(fit)
abline(fit8, col=2, lty=2)
abline(fit16, col=3, lty=2)
coef8 <- coefficients(fit8)
coef16 <- coefficients(fit16)
sig <- 3
text(48,30,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=.8)
text(48,27,labels = paste("16 mm: ",signif(coef16[2],sig),"x+",signif(coef16[1],sig), sep=" "),col=3, cex=.8)
legend("topleft",legend = c("Proportionerede styrker","Max. tilslag: 8 mm", "Max. tilslag: 16 mm"), lty=c(5,0,0), pch=c(NA,16,16), lwd=1.5, col=c(1,2,3))
# ------------------------------------------------------------------
par(pty="s")
plot(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=2)
grid(equilogs = T, col="darkgrey")
abline(v=prop, lty=5, lwd=1.5)
Plot_ErrorBars(x = mutryk[,1],y = mulokc[,1],x.err = sdtryk[,1],y.err = sdlokc[,1],xbar = T,ybar = T)
Plot_ErrorBars(x = mutryk[,2],y = mulokc[,2],x.err = sdtryk[,2],y.err = sdlokc[,2],xbar = T,ybar = T)
points(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,40), ylim=c(10,40), pch=16, col=2)
points(mutryk[,2],mulokc[,2], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=3)
fit8 <- lm(mulokc[,1]~mutryk[,1])
fit16 <- lm(mulokc[,2]~mutryk[,2])
summary(fit)
abline(fit8, col=2, lty=2)
abline(fit16, col=3, lty=2)
coef8 <- coefficients(fit8)
coef16 <- coefficients(fit16)
sig <- 3
text(48,30,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=.8)
text(48,27,labels = paste("16 mm: ",signif(coef16[2],sig),"x+",signif(coef16[1],sig), sep=" "),col=3, cex=.8)
legend(locator(1),legend = c("Proportionerede styrker","Max. tilslag: 8 mm", "Max. tilslag: 16 mm"), lty=c(5,0,0), pch=c(NA,16,16), lwd=1.5, col=c(1,2,3))
par(pty="s",xpd=F)
plot(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=2)
grid(equilogs = T, col="darkgrey")
abline(v=prop, lty=5, lwd=1.5)
Plot_ErrorBars(x = mutryk[,1],y = mulokc[,1],x.err = sdtryk[,1],y.err = sdlokc[,1],xbar = T,ybar = T)
Plot_ErrorBars(x = mutryk[,2],y = mulokc[,2],x.err = sdtryk[,2],y.err = sdlokc[,2],xbar = T,ybar = T)
points(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,40), ylim=c(10,40), pch=16, col=2)
points(mutryk[,2],mulokc[,2], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=3)
fit8 <- lm(mulokc[,1]~mutryk[,1])
fit16 <- lm(mulokc[,2]~mutryk[,2])
summary(fit)
abline(fit8, col=2, lty=2)
abline(fit16, col=3, lty=2)
coef8 <- coefficients(fit8)
coef16 <- coefficients(fit16)
sig <- 3
text(48,30,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=.8)
text(48,27,labels = paste("16 mm: ",signif(coef16[2],sig),"x+",signif(coef16[1],sig), sep=" "),col=3, cex=.8)
legend(locator(1),legend = c("Proportionerede styrker","Max. tilslag: 8 mm", "Max. tilslag: 16 mm"), lty=c(5,0,0), pch=c(NA,16,16), lwd=1.5, col=c(1,2,3))
par(pty="s",xpd=F)
plot(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=2)
grid(equilogs = T, col="darkgrey")
abline(v=prop, lty=5, lwd=1.5)
Plot_ErrorBars(x = mutryk[,1],y = mulokc[,1],x.err = sdtryk[,1],y.err = sdlokc[,1],xbar = T,ybar = T)
Plot_ErrorBars(x = mutryk[,2],y = mulokc[,2],x.err = sdtryk[,2],y.err = sdlokc[,2],xbar = T,ybar = T)
points(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,40), ylim=c(10,40), pch=16, col=2)
points(mutryk[,2],mulokc[,2], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=3)
fit8 <- lm(mulokc[,1]~mutryk[,1])
fit16 <- lm(mulokc[,2]~mutryk[,2])
summary(fit)
abline(fit8, col=2, lty=2)
abline(fit16, col=3, lty=2)
coef8 <- coefficients(fit8)
coef16 <- coefficients(fit16)
sig <- 3
text(48,30,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=.8)
text(48,27,labels = paste("16 mm: ",signif(coef16[2],sig),"x+",signif(coef16[1],sig), sep=" "),col=3, cex=.8)
legend("topright"",legend = c("Proportionerede styrker","Max. tilslag: 8 mm", "Max. tilslag: 16 mm"), lty=c(5,0,0), pch=c(NA,16,16), lwd=1.5, col=c(1,2,3))
# ----
legend("topright",legend = c("Proportionerede styrker","Max. tilslag: 8 mm", "Max. tilslag: 16 mm"), lty=c(5,0,0), pch=c(NA,16,16), lwd=1.5, col=c(1,2,3))
par(pty="s",xpd=F)
plot(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=2)
grid(equilogs = T, col="darkgrey")
abline(v=prop, lty=5, lwd=1.5)
Plot_ErrorBars(x = mutryk[,1],y = mulokc[,1],x.err = sdtryk[,1],y.err = sdlokc[,1],xbar = T,ybar = T)
Plot_ErrorBars(x = mutryk[,2],y = mulokc[,2],x.err = sdtryk[,2],y.err = sdlokc[,2],xbar = T,ybar = T)
points(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,40), ylim=c(10,40), pch=16, col=2)
points(mutryk[,2],mulokc[,2], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=3)
fit8 <- lm(mulokc[,1]~mutryk[,1])
fit16 <- lm(mulokc[,2]~mutryk[,2])
summary(fit)
abline(fit8, col=2, lty=2)
abline(fit16, col=3, lty=2)
coef8 <- coefficients(fit8)
coef16 <- coefficients(fit16)
sig <- 3
text(48,30,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=.8)
text(48,27,labels = paste("16 mm: ",signif(coef16[2],sig),"x+",signif(coef16[1],sig), sep=" "),col=3, cex=.8)
legend("topleft",legend = c("Proportionerede styrker","Max. tilslag: 8 mm", "Max. tilslag: 16 mm"), lty=c(5,0,0), pch=c(NA,16,16), lwd=1.5, col=c(1,2,3))
plot(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=2)
grid(equilogs = T, col="darkgrey")
abline(v=prop, lty=5, lwd=1.5)
Plot_ErrorBars(x = mutryk[,1],y = mulokc[,1],x.err = sdtryk[,1],y.err = sdlokc[,1],xbar = T,ybar = T)
Plot_ErrorBars(x = mutryk[,2],y = mulokc[,2],x.err = sdtryk[,2],y.err = sdlokc[,2],xbar = T,ybar = T)
points(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,40), ylim=c(10,40), pch=16, col=2)
points(mutryk[,2],mulokc[,2], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=3)
fit8 <- lm(mulokc[,1]~mutryk[,1])
fit16 <- lm(mulokc[,2]~mutryk[,2])
summary(fit)
abline(fit8, col=2, lty=2)
abline(fit16, col=3, lty=2)
coef8 <- coefficients(fit8)
coef16 <- coefficients(fit16)
sig <- 3
text(48,30,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=.8)
text(48,27,labels = paste("16 mm: ",signif(coef16[2],sig),"x+",signif(coef16[1],sig), sep=" "),col=3, cex=.8)
legend("topleft",legend = c("Proportionerede styrker","Max. tilslag: 8 mm", "Max. tilslag: 16 mm"), lty=c(5,0,0), pch=c(NA,16,16), lwd=1.5, col=c(1,2,3))
# ------------------
beton
View(beton)
mulokc
sdlokc
beton
beton[beton$prop.styrke==34,]
length(beton[beton$prop.styrke==34,])
str(beton[beton$prop.styrke==34,])
length(beton[beton$prop.styrke==20,])
beton[beton$prop.styrke==20,]
beton[beton$prop.styrke==27,]
beton[beton$prop.styrke==40,]
beton[beton$prop.styrke==40,]
signif(cbind(mutryk[,2],sdtryk[,1],mulokc[,2],sdlokc[,1]),3)
red <- edata[-c(1:3),]
plot(sigma~tojning, data=edata, xlab="toejning [%]", ylab="Spaending [MPa]", xlim=c(0,max(edata$tojning)), ylim=c(0,max(edata$sigma)), cex=1.2)
fit <- lm(sigma~tojning, data=edata)
summary(fit)
abline(fit)
points(sigma~tojning, data=red, pch=16, col=2, cex=0.9)
fit2 <- lm(sigma~tojning, data=red)
summary(fit2)
abline(fit2,col=2)
sig <- 5
text(x = 5e-3,y = c(15,13),labels = c(paste("E: ", signif(coefficients(fit)[2]*100,sig), " MPa", sep=""),paste("E: ", signif(coefficients(fit2)[2]*100,sig), " MPa", sep="")),col=c(1,2))
text(x = c(5e-3,3e-2),y = c(0.4,0.1)*brudcyl+1,labels = c("40% af brud","10% af brud"), col="darkgrey")
abline(h=c(0.1*brudcyl,0.4*brudcyl), lty=2, col="darkgrey")
text(x = 5e-3,y = 13,labels = paste("E: ", signif(coefficients(fit2)[2]*100,sig), " MPa", sep=""),col=2)
plot(sigma~tojning, data=edata, xlab="toejning [%]", ylab="Spaending [MPa]", xlim=c(0,max(edata$tojning)), ylim=c(0,max(edata$sigma)), cex=1.2)
fit <- lm(sigma~tojning, data=edata)
summary(fit)
abline(fit)
points(sigma~tojning, data=red, pch=16, col=2, cex=0.9)
fit2 <- lm(sigma~tojning, data=red)
summary(fit2)
abline(fit2,col=2)
sig <- 5
text(x = 5e-3,y = 13,labels = paste("E: ", signif(coefficients(fit2)[2]*100,sig), " MPa", sep=""),col=2)
text(x = c(5e-3,3e-2),y = c(0.4,0.1)*brudcyl+1,labels = c("40% af brud","10% af brud"), col="darkgrey")
abline(h=c(0.1*brudcyl,0.4*brudcyl), lty=2, col="darkgrey")
plot(sigma~tojning, data=edata, xlab="toejning [%]", ylab="Spaending [MPa]", xlim=c(0,max(edata$tojning)), ylim=c(0,max(edata$sigma)), cex=1.2)
fit <- lm(sigma~tojning, data=edata)
summary(fit)
abline(fit)
points(sigma~tojning, data=red, pch=16, col=2, cex=0.9)
fit2 <- lm(sigma~tojning, data=red)
summary(fit2)
abline(fit2,col=2)
sig <- 3
text(x = 5e-3,y = 13,labels = paste("E: ", signif(coefficients(fit2)[2]*100,sig), " MPa", sep=""),col=2)
text(x = c(5e-3,3e-2),y = c(0.4,0.1)*brudcyl+1,labels = c("40% af brud","10% af brud"), col="darkgrey")
abline(h=c(0.1*brudcyl,0.4*brudcyl), lty=2, col="darkgrey")
coefficients(fit2)[2]
text(x = 5e-3,y = 13,labels = paste("E: ", signif(coefficients(fit2)[2]/10,sig), " MPa", sep=""),col=2)
plot(sigma~tojning, data=edata, xlab="toejning [%]", ylab="Spaending [MPa]", xlim=c(0,max(edata$tojning)), ylim=c(0,max(edata$sigma)), cex=1.2)
fit <- lm(sigma~tojning, data=edata)
summary(fit)
abline(fit)
points(sigma~tojning, data=red, pch=16, col=2, cex=0.9)
fit2 <- lm(sigma~tojning, data=red)
summary(fit2)
abline(fit2,col=2)
sig <- 3
text(x = 5e-3,y = 13,labels = paste("E: ", signif(coefficients(fit2)[2]/10,sig), " GPa", sep=""),col=2)
text(x = c(5e-3,3e-2),y = c(0.4,0.1)*brudcyl+1,labels = c("40% af brud","10% af brud"), col="darkgrey")
abline(h=c(0.1*brudcyl,0.4*brudcyl), lty=2, col="darkgrey")
beton8
par(pty="s",xpd=F)
plot(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=2)
grid(equilogs = T, col="darkgrey")
abline(v=prop, lty=5, lwd=1.5)
Plot_ErrorBars(x = mutryk[,1],y = mulokc[,1],x.err = sdtryk[,1],y.err = sdlokc[,1],xbar = T,ybar = T)
Plot_ErrorBars(x = mutryk[,2],y = mulokc[,2],x.err = sdtryk[,2],y.err = sdlokc[,2],xbar = T,ybar = T)
points(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,40), ylim=c(10,40), pch=16, col=2)
points(mutryk[,2],mulokc[,2], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=3)
fit8 <- lm(mulokc[,1]~mutryk[,1])
fit16 <- lm(mulokc[,2]~mutryk[,2])
summary(fit)
abline(fit8, col=2, lty=2)
abline(fit16, col=3, lty=2)
coef8 <- coefficients(fit8)
coef16 <- coefficients(fit16)
sig <- 3
text(48,30,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=.8)
text(48,27,labels = paste("16 mm: ",signif(coef16[2],sig),"x+",signif(coef16[1],sig), sep=" "),col=3, cex=.8)
legend("topleft",legend = c("Proportionerede styrker","Max. tilslag: 8 mm", "Max. tilslag: 16 mm"), lty=c(5,0,0), pch=c(NA,16,16), lwd=1.5, col=c(1,2,3))
# ---
text(48,30,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=1.1)
text(48,27,labels = paste("16 mm: ",signif(coef16[2],sig),"x+",signif(coef16[1],sig), sep=" "),col=3, cex=1.1)
text(45,30,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=1.1)
par(pty="s",xpd=F)
plot(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=2)
grid(equilogs = T, col="darkgrey")
abline(v=prop, lty=5, lwd=1.5)
Plot_ErrorBars(x = mutryk[,1],y = mulokc[,1],x.err = sdtryk[,1],y.err = sdlokc[,1],xbar = T,ybar = T)
Plot_ErrorBars(x = mutryk[,2],y = mulokc[,2],x.err = sdtryk[,2],y.err = sdlokc[,2],xbar = T,ybar = T)
points(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,40), ylim=c(10,40), pch=16, col=2)
points(mutryk[,2],mulokc[,2], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=3)
fit8 <- lm(mulokc[,1]~mutryk[,1])
fit16 <- lm(mulokc[,2]~mutryk[,2])
summary(fit)
abline(fit8, col=2, lty=2)
abline(fit16, col=3, lty=2)
coef8 <- coefficients(fit8)
coef16 <- coefficients(fit16)
sig <- 3
text(45,30,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=1.1)
text(45,27,labels = paste("16 mm: ",signif(coef16[2],sig),"x+",signif(coef16[1],sig), sep=" "),col=3, cex=1.1)
legend("topleft",legend = c("Proportionerede styrker","Max. tilslag: 8 mm", "Max. tilslag: 16 mm"), lty=c(5,0,0), pch=c(NA,16,16), lwd=1.5, col=c(1,2,3))
# -------
text(45,25,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=1.1)
text(45,23,labels = paste("16 mm: ",signif(coef16[2],sig),"x+",signif(coef16[1],sig), sep=" "),col=3, cex=1.1)
par(pty="s",xpd=F)
plot(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=2)
grid(equilogs = T, col="darkgrey")
abline(v=prop, lty=5, lwd=1.5)
Plot_ErrorBars(x = mutryk[,1],y = mulokc[,1],x.err = sdtryk[,1],y.err = sdlokc[,1],xbar = T,ybar = T)
Plot_ErrorBars(x = mutryk[,2],y = mulokc[,2],x.err = sdtryk[,2],y.err = sdlokc[,2],xbar = T,ybar = T)
points(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,40), ylim=c(10,40), pch=16, col=2)
points(mutryk[,2],mulokc[,2], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=3)
fit8 <- lm(mulokc[,1]~mutryk[,1])
fit16 <- lm(mulokc[,2]~mutryk[,2])
summary(fit)
abline(fit8, col=2, lty=2)
abline(fit16, col=3, lty=2)
coef8 <- coefficients(fit8)
coef16 <- coefficients(fit16)
sig <- 3
text(45,25,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=1.1)
text(45,22,labels = paste("16 mm: ",signif(coef16[2],sig),"x+",signif(coef16[1],sig), sep=" "),col=3, cex=1.1)
legend("topleft",legend = c("Proportionerede styrker","Max. tilslag: 8 mm", "Max. tilslag: 16 mm"), lty=c(5,0,0), pch=c(NA,16,16), lwd=1.5, col=c(1,2,3))
# ----
par(pty="s",xpd=F)
plot(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=2)
grid(equilogs = T, col="darkgrey")
abline(v=prop, lty=5, lwd=1.5)
Plot_ErrorBars(x = mutryk[,1],y = mulokc[,1],x.err = sdtryk[,1],y.err = sdlokc[,1],xbar = T,ybar = T)
Plot_ErrorBars(x = mutryk[,2],y = mulokc[,2],x.err = sdtryk[,2],y.err = sdlokc[,2],xbar = T,ybar = T)
points(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,40), ylim=c(10,40), pch=16, col=2)
points(mutryk[,2],mulokc[,2], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=3)
fit8 <- lm(mulokc[,1]~mutryk[,1])
fit16 <- lm(mulokc[,2]~mutryk[,2])
summary(fit)
abline(fit8, col=2, lty=2)
abline(fit16, col=3, lty=2)
coef8 <- coefficients(fit8)
coef16 <- coefficients(fit16)
sig <- 3
text(45,25,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=1.1)
text(45,22,labels = paste("16 mm: ",signif(coef16[2],sig),"x+",signif(coef16[1],sig), sep=" "),col=3, cex=1.1)
legend("topleft",legend = c("Proportionerede styrker","Max. tilslag: 8 mm", "Max. tilslag: 16 mm"), lty=c(5,0,0), pch=c(NA,16,16), lwd=1.5, col=c(1,2,3))
# ----------------
text(40,20,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=1.1)
text(40,17,labels = paste("16 mm: ",signif(coef16[2],sig),"x+",signif(coef16[1],sig), sep=" "),col=3, cex=1.1)
par(pty="s",xpd=F)
plot(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=2)
grid(equilogs = T, col="darkgrey")
abline(v=prop, lty=5, lwd=1.5)
Plot_ErrorBars(x = mutryk[,1],y = mulokc[,1],x.err = sdtryk[,1],y.err = sdlokc[,1],xbar = T,ybar = T)
Plot_ErrorBars(x = mutryk[,2],y = mulokc[,2],x.err = sdtryk[,2],y.err = sdlokc[,2],xbar = T,ybar = T)
points(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,40), ylim=c(10,40), pch=16, col=2)
points(mutryk[,2],mulokc[,2], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=3)
fit8 <- lm(mulokc[,1]~mutryk[,1])
fit16 <- lm(mulokc[,2]~mutryk[,2])
summary(fit)
abline(fit8, col=2, lty=2)
abline(fit16, col=3, lty=2)
coef8 <- coefficients(fit8)
coef16 <- coefficients(fit16)
sig <- 3
text(40,20,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=1.1)
text(40,17,labels = paste("16 mm: ",signif(coef16[2],sig),"x+",signif(coef16[1],sig), sep=" "),col=3, cex=1.1)
legend("topleft",legend = c("Proportionerede styrker","Max. tilslag: 8 mm", "Max. tilslag: 16 mm"), lty=c(5,0,0), pch=c(NA,16,16), lwd=1.5, col=c(1,2,3))
# -------------------
plot(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=2)
grid(equilogs = T, col="darkgrey")
abline(v=prop, lty=5, lwd=1.5)
Plot_ErrorBars(x = mutryk[,1],y = mulokc[,1],x.err = sdtryk[,1],y.err = sdlokc[,1],xbar = T,ybar = T)
Plot_ErrorBars(x = mutryk[,2],y = mulokc[,2],x.err = sdtryk[,2],y.err = sdlokc[,2],xbar = T,ybar = T)
points(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,40), ylim=c(10,40), pch=16, col=2)
points(mutryk[,2],mulokc[,2], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=3)
fit8 <- lm(mulokc[,1]~mutryk[,1])
fit16 <- lm(mulokc[,2]~mutryk[,2])
summary(fit)
abline(fit8, col=2, lty=2)
abline(fit16, col=3, lty=2)
coef8 <- coefficients(fit8)
coef16 <- coefficients(fit16)
sig <- 3
text(40,20,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=1.1)
text(40,17,labels = paste("16 mm: ",signif(coef16[2],sig),"x+",signif(coef16[1],sig), sep=" "),col=3, cex=1.1)
legend("topleft",legend = c("Proportionerede styrker","Max. tilslag: 8 mm", "Max. tilslag: 16 mm"), lty=c(5,0,0), pch=c(NA,16,16), lwd=1.5, col=c(1,2,3))
#
plot(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=2)
grid(equilogs = T, col="darkgrey")
abline(v=prop, lty=5, lwd=1.5)
Plot_ErrorBars(x = mutryk[,1],y = mulokc[,1],x.err = sdtryk[,1],y.err = sdlokc[,1],xbar = T,ybar = T)
Plot_ErrorBars(x = mutryk[,2],y = mulokc[,2],x.err = sdtryk[,2],y.err = sdlokc[,2],xbar = T,ybar = T)
points(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,40), ylim=c(10,40), pch=16, col=2)
points(mutryk[,2],mulokc[,2], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=3)
fit8 <- lm(mulokc[,1]~mutryk[,1])
fit16 <- lm(mulokc[,2]~mutryk[,2])
summary(fit)
abline(fit8, col=2, lty=2)
abline(fit16, col=3, lty=2)
coef8 <- coefficients(fit8)
coef16 <- coefficients(fit16)
sig <- 3
text(40,20,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=1.1)
text(40,17,labels = paste("16 mm: ",signif(coef16[2],sig),"x+",signif(coef16[1],sig), sep=" "),col=3, cex=1.1)
legend("topleft",legend = c("Proportionerede styrker","Max. tilslag: 8 mm", "Max. tilslag: 16 mm"), lty=c(5,0,0), pch=c(NA,16,16), lwd=1.5, col=c(1,2,3))
#
text(41,20,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=1.1)
plot(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=2)
grid(equilogs = T, col="darkgrey")
abline(v=prop, lty=5, lwd=1.5)
Plot_ErrorBars(x = mutryk[,1],y = mulokc[,1],x.err = sdtryk[,1],y.err = sdlokc[,1],xbar = T,ybar = T)
Plot_ErrorBars(x = mutryk[,2],y = mulokc[,2],x.err = sdtryk[,2],y.err = sdlokc[,2],xbar = T,ybar = T)
points(mutryk[,1],mulokc[,1], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,40), ylim=c(10,40), pch=16, col=2)
points(mutryk[,2],mulokc[,2], xlab="Trykstyrke MPa", ylab="Lokstyrke MPa", xlim=c(10,max(mulokc)+5), ylim=c(10,max(mutryk)+5), pch=16, col=3)
fit8 <- lm(mulokc[,1]~mutryk[,1])
fit16 <- lm(mulokc[,2]~mutryk[,2])
summary(fit)
abline(fit8, col=2, lty=2)
abline(fit16, col=3, lty=2)
coef8 <- coefficients(fit8)
coef16 <- coefficients(fit16)
sig <- 3
text(41,20,labels = paste("8 mm:   ",signif(coef8[2],sig),"x+",signif(coef8[1],sig), sep=" "),col=2, cex=1.1)
text(41,17,labels = paste("16 mm: ",signif(coef16[2],sig),"x+",signif(coef16[1],sig), sep=" "),col=3, cex=1.1)
legend("topleft",legend = c("Proportionerede styrker","Max. tilslag: 8 mm", "Max. tilslag: 16 mm"), lty=c(5,0,0), pch=c(NA,16,16), lwd=1.5, col=c(1,2,3))
2+2
hyper
phyper(q = 4,m = 20,n = 30, k=10)
210*20^4*30^6/(50^10)
dhyper(q = 4,m = 20,n = 30, k=10)
phyper(q = 4,m = 20,n = 30, k=10)
phyper(q = 4,m = 20,n = 30, k=10)-phyper(q = 3,m = 20,n = 30, k=10)
dhyper(x=4,m=20,n=30,k=10)
10!
factorial(3)
factorial(20)/(factorial(4)*factorial(16))*factorial(30)(factorial(6)*factorial(24))/(factorial(50)/(factorial(10)*factorial(40)))
factorial(20)/(factorial(4)*factorial(16))*factorial(30)(factorial(6)*factorial(24))/(factorial(50)/(factorial(10)*factorial(40)))
(factorial(20)/(factorial(4)*factorial(16))*factorial(30)/(factorial(6)*factorial(24)))/(factorial(50)/(factorial(10)*factorial(40)))
phyper(q = 4,m = 20,n = 30,k = 10)
phyper(q = 4,m = 20,n = 30,k = 10)-phyper(q = 3,m = 20,n = 30,k = 10)
pbinom(q = 4,size = 10,prob = 20/50)
210*20^4*30^6/50^10
pbinom(q = 4,size = 10,prob = 20/50)-pbinom(q = 3,size = 10,prob = 20/50)
dbinom(x = 4,size = 10,prob = 20/50)
dhyper(x = 4,m = 20,n = 30,k=10)
ph <- c(0.5, 0.45, 0.5)
ph
sum(ph)
ph <- c(0.5, 0.45, 0.05)
ph
sum(ph)
ph <- c(0.5, 0.45, 0.05)
(ph <- c(0.5, 0.45, 0.05))
(ph <- c(0.5, 0.45, 0.05))
pah <- c(0.10, 0.01, 0.39)
PHA <- c()
PHA(1)
PHA(1) <- 1
PHA[1] <- 1
PHA
PHA[2] <- 3
PHA
PHA <- c()
PHA
bund <- 0
bund <- 0
for(i in 1:3){
bund <- bund + ph[i]*pah[i]
}
bund
for(i in 1:3){
PHA[i] <- pah[i]*ph[i]/bund
}
PHA
sum(PHA)
dbinom
pbinom(q = 5000,size = 10000, prob = 0.5)
dbinom(q = 5000,size = 10000, prob = 0.5)
dbinom(x = 5000,size = 10000, prob = 0.5)
qbinom(p = 1/6,size = 10000, prob = 0.5)
qbinom(p = 5/6,size = 10000, prob = 0.5)
qbinom(p = 1/2,size = 10000, prob = 0.5)
qbinom(p = 5/6,size = 10000, prob = 0.5)
qbinom(p = 0.475,size = 10000, prob = 0.5)
qbinom(p = 0.025,size = 10000, prob = 0.5)
qbinom(p = 0.975,size = 10000, prob = 0.5)
qbinom(p = 1/6,size = 10000, prob = 0.5)
qbinom(p = 5/6,size = 10000, prob = 0.5)
sqrt(10000*0.5*0.5)
qbinom(p = 1/6,size = 0, prob = 0.5)
qbinom(p = 1/6,size = 100, prob = 0.5)
qbinom(p = 1/6,size = 10000, prob = 0.5)
qbinom(p = 1/6,size = 100, prob = 0.5)
qbinom(p = 1/6,size = 1000, prob = 0.5)
qbinom(p = 1/6,size = 10000, prob = 0.5)
sd(52,60,55,52,54,52,46,42,50,54,54,54)
sd(c(52,60,55,52,54,52,46,42,50,54,54,54))
dnorm(x = -0,19)
pnorm(q = 9.8, mean = 9.78, sd = 0.0031)
pnorm(q = 9.784, mean = 9.78, sd = 0.0031)
1-pnorm(q = 9.784, mean = 9.78, sd = 0.0031)
pnorm(q = 9.7794, mean = 9.78, sd = 0.0031)
qnorm(p = 0.9,mean = 9.78, sd = 0.0031)
qnorm(p = 0.91,mean = 9.78, sd = 0.0031, )
qnorm(p = 0.9,mean = 9.78, sd = 0.0031)
pnorm(q = 9.8, mean = 9.78, sd = 0.0031) - pnorm(q = 9.784, mean = 9.78, sd = 0.0031)
pnorm(q = 9.7794, mean = 9.78, sd = 0.0031)
setwd("~/Google Drive/Kangaamiut/Feltundersogelser/UCS/Data logger")
(folder <- list.dirs())
(Pname <- list.files())
pmax
Pmax
(folder <- list.dirs())
(Pname <- list.files())
# Konvertering fra Volt til kN
convfact <- 50 # kN/V
library("lubridate")
Pmax <- c()
tlim <- data.frame()
# Bestem Pmax og tlim
plim <- c(NA,rep(10,23))
plim[6] <- 12
plim[20:24] <- 45
for(i in 2:24){
input = folder[i]
file <-paste(input,"/",list.files(input)[1],sep = "")
(txt <- read.csv(file, sep = ",", header = FALSE, skip = 27))
# Kolonnenavne
colnames(txt) <- c("Nr", "Time", "VDC", "Limit")
# Index
#index = which(folder == input)
# Last [kN]
txt$Load <- txt$VDC*convfact
# Tid
t <- as.POSIXct(txt$Time, format = "%d-%m-%Y %H:%M:%OS")
txt$Tid <- minute(t)*60 + second(t) - minute(t[1])*60
# Max last [kN]
pmax <- max(txt$Load)
Pmax[i] <- pmax
tmax <- txt$Tid[which(txt$Load == max(txt$Load))][1]
tlim[i,1] <- txt$Tid[min(which(txt$Load > plim[i]))]
tlim[i,2] <- tmax
}
txt
par(mfrow = c(5,5), mar = c(2,2,1,1))
#myplot <- "pmax_graf"
#pdf(paste0("~/6931816wdyfbsdksggh/", myplot, ".pdf"))
LoadR <- c()
for(i in 2:24){
input = folder[i]
file <-paste(input,"/",list.files(input)[1],sep = "")
(txt <- read.csv(file, sep = ",", header = FALSE, skip = 27))
# Kolonnenavne
colnames(txt) <- c("Nr", "Time", "VDC", "Limit")
# Last [kN]
txt$Load <- txt$VDC*convfact
# Tid
t <- as.POSIXct(txt$Time, format = "%d-%m-%Y %H:%M:%OS")
txt$Tid <- minute(t)*60 + second(t) - minute(t[1])*60
# Plot
plot(Load~Tid, txt, type = "l")
points(tlim[i,2],Pmax[i], col = 2, pch = 16)
title(main = paste(folder[i]))
text(x = tlim[i,2] - 15, y = Pmax[i]-5, labels = round(Pmax[i],1))
abline(v = tlim[i,1], col = 3)
abline(v = tlim[i,2], col = 2)
Lrate <- lm(Load~Tid,txt[tlim[i,1] < txt$Tid & txt$Tid < tlim[i,2],])
LoadR[i] <- coefficients(Lrate)[2]
abline(Lrate, col = 4)
}
